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Abstract 

We calculate the rapidity dependence of J /if/ nuclear modification factor and averaged transverse 
momentum square in heavy ion collisions at RHIC in a 3-dimensional transport approach with 
regeneration mechanism. 
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J I 'if/ suppression flJL£| is widely accepted as a probe of quark-gluon plasma (QGP) formed in 
relativistic heavy ion collisions and was first observed at SPS [2] more than ten years ago. At 
RHIC and LHC energy, a significant number of charm quarks are generated in central heavy ion 
collisions, and the recombination of these uncorrelated charm quarks offers another source for 
J /if/ production [3]. There are different ways to describe the J/i/i regeneration. In the statisti- 
cal model |4|], all the J/if/s are produced at hadronization of the system through thermal distri- 
butions and charm conservation. In some other models, J/i//s come from both the continuous 
regeneration inside the hot medium and primordial production through initial nucleon-nucleon 
collisions 

ma. 

The regeneration is used to describe [6] the J /if/ nuclear modification factor 
Raa and averaged transverse momentum square (p 2 ). From the experimental data |7[ at RHIC, 
the Raa is almost the same as that at SPS, which seems difficult to explain in models with only 
primordial production mechanism. 

The rapidity dependence of J/ip production was also measured at RHIC 10. Si and discussed 
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in models 12, [KJ, 11 1H . The surprising finding of the experimental result is that the apparent 
suppression at forward rapidity is stronger than that at midrapidity, i.e. > R^ Hard . This is 
again difficult to explain in models with only initial production mechanism, since the suppression 
at midrapidity should be stronger than that at forward rapidity. Not only Raa but also (/? 2 ) 
depends on the rapidity |7|]. In semi-central and central Au+Au collisions the value of (pf) at 
midrapidity is lower than that at forward rapidity, i.e. {p\)' md < {p^)f orward . In this paper, we 
will discuss the rapidity dependence of Raa and (p 2 ) in a 3-dimensional transport model with 
both initial production and continuous regeneration mechanisms. 

At RHIC J/if/s are detected at midrapidity \y\ < 0.35 and forward rapidity 1.2 < y < 2.2, 
both are located in the plateau of the rapidity distribution of light hadrons B12I1 . Therefore, the 
space-time evolution of the QGP can be approximately described by the transverse hydrodynamic 
equations at midrapidity 11311 . and the J/if/ motion is controlled by a 3-dimensional transport 
equation in an explicitly boost invariant form 

[cosh^ - 77) d/dr + sinhfjy - t])/t d/di] + ■ V,J = -a^ fa + (3^, (1) 
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where fa = fy(Pi,y, x t , rj, r|b) is the charmonium distribution function in phase space at fixed im- 
pact parameter b, and we have used transverse energy E t = ^E 2 ^, — p 2 , rapidity = 1/2 ln[(ity+ 
p z )/(Ey - p z )], proper time t = Vf 2 - z 2 and space-time rapidity rj — 1/2 In [(f + z)/(f - z)] to 
replace the charmonium energy £V = -^p 2 + m 2 , longitudinal momentum p,, time t and lon- 
gitudinal coordinate z- The term with transverse velocity = p,/E, reflects the leakage effect 
in charmonium motion. To take into account the decay of the charmonium excitation states into 
J/if/, the symbol *P here stands for J/tf/,x c and iff' and the ratio of their contributions in the ini- 
tial condition is taken as 6:3:1. The suppression and regeneration in the QGP are described by 
the loss and gain terms on the right hand side of the transport equation. Considering the gluon 



dissociation process, a can be explicitly written as 111411 



Mp, x, f|b) = J d\l ({2n?AE,E g ) W%{s)f g {q, T, u)&(T - T c )/®(Tj - T), (2) 

where W ^ is the transition probability lfl5ll as a function of the colliding energy yfs of the disso- 
ciation process, E g and f g are the gluon energy and gluon thermal distribution, and T c and rj are 
the critical temperature of the deconfinement phase transition and dissociation temperature of *P, 
taken as T c = 165 MeV, T ] J* /T c = 1.9 and Tf /T c = T$ /T c = 1. The two step functions in 
a indicate that the suppression is finite in the QGP phase at temperature T < Tj and becomes 
infinite at T > Tj, We have here neglected the suppression process in hadron phase lfl3i l6ll. 

The gain term /3y can be obtained from the loss term a by considering detailed balance fl. 
We assume local thermalization of charm quarks in the QGP and take the charm quark distribu- 
tion as 

f c (k,x,t)=p c (x,t)f q (k) (3) 
with p c being the density of charm quarks in coordinate space, 

p c (x, t) = T A (x,)T B (x, - b) cosh tj/t da^ N ldr\ (4) 

and f q the normalized Fermi distribution in momentum space, where 7^ and Tb are the thickness 
functions for the two colliding nuclei determined by nuclear geometry. Since the large uncer- 
tainty of charm quark production cross section in pp collisions for both experimental and theo- 
retical studies, we assume the rapidity dependence of charm production as a Gauss distribution 
do- C p p /dr] = da^ p /dT}\ =Q e~'' l2l %> with do^p/dr]^^ = 120 pb which agrees with the experimental 

data iflill and (dcr^pldr^ ^ ^/(dcrp'p/di]^^) = 1/3 to determine the parameter tjq which is in 
between the smallest and largest theoretical estimation lfl6Tl . 

The contribution from the primordial charmonium production is reflected in the initial con- 
dition of the transport equation at the starting time to. By fitting the experimental data ifnll for 
pp collisions at RHIC, the initial charmonium momentum distribution is extracted as 

f pp {Vny) = 5g(y)/ {An{p 2 )(y)) [l + p 2 / (4<^>(y))f 6 , (5) 



where the rapidity distribution g{y) is a double Gauss function 1 17], and the rapidity dependence 
of the averaged transverse momentum square is taken as (pf)(y) - (p 2 )(0)(l - y 2 ly 2 wi ) with the 
parameters (pf)(0) =4.1 (GeV/c) 2 and y mflv = arccosh( yfs / 2m j /,/,). Note that, to include the 
Cronin effect in the initial state of heavy ion collisions 11811. we add an extra term to (p 2 ) which 
comes from the gluon multi-scattering with nucleons 1161 1 1 ill . 
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The charmonium production, including initial production and regeneration, is related to the 
QGP evolution through the local temperature T and fluid velocity appearing in the thermal 
gluon and charm quark distributions, they are determined by the ideal hydrodynamics II 1311 . 
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Figure 1: The nuclear modification factor Raa (left panel) and averaged transverse momentum square (pj) (right panel) 
at mid and forward rapidity as functions of number of participants N p . The theoretical calculations with only initial 
production (dot-dashed lines), only regeneration (dashed lines) and both (solid lines) are compared with the experimental 
data 01. 



With the known distribution fj/,/,(p t ,y,Xf, 7], r|b), one can calculate the J /if/ yield and momen- 
tum spectra. The nuclear modification factor Raa and averaged transverse momentum square (pf) 
at mid and forward rapidity are shown in FigQ]as functions of centrality. Since Raa is normal- 
ized to the pp collisions, the assumption of the same medium at mid and forward rapidity leads 
to similar Raa in the two rapidity regions, when we consider only initial production, as shown 
in the left panel. The regeneration at forward rapidity is, however, much less than that at mid 
rapidity. As a result of the competition, the total Raa at forward rapidity is less than that at mid 
rapidity, consistent with the experimental observation. 

While the population is dominated by low momentum J/if/s, the averaged transverse momen- 
tum carries more information on high momentum J/if/s and can tell us more about the dynam- 
ics of charmonium production and suppression. The initially produced J/if/s are from the hard 
nucleon-nucleon process at the very beginning of the collision and their p, spectrum is harder. 
From the gluon multi-scattering with nucleons before the two gluons fuse into a J /if/, there is 
a p t broadening for the initially produced J/i//s. Considering further the leakage effect which 
enables the high momentum J/iJ/s, escape from the anomalous suppression in the hot medium, 
the initially produced (pf) increases smoothly with centrality and becomes saturated at large N p . 
Since the regenerated J/if/s are from the thermalized charm quarks inside the QGP, their averaged 
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momentum is small and almost independent of the centrality. Both the initially produced and re- 
generated (pf) is not sensitive to the rapidity region. While the difference between the initially 
produced and regenerated Raa decreases with increasing N p , the difference between the values 
of (pf) from the two rapidity regions increases smoothly with centrality! The total (pf) depends 
strongly on the fraction of the regeneration. At mid rapidity, the regeneration and initial produc- 
tion are equally important in central collisions, see the left panel of Fig[T] The large contribution 
from the regeneration leads to a remarkable decrease of the value of (pf) at mid rapidity. At for- 
ward rapidity, the regeneration contribution is, however, very small even for central collisions, 
see the left panel again. In this case, the total (pf) is dominated by the initial production in the 
whole N p region. 

In summary, we calculated the J/iJ/ nuclear modification factor and averaged transverse mo- 
mentum square at mid and forward rapidity in a three dimensional transport approach. The 
experimentally observed rapidity dependence of Raa and (pf) in Au+Au collisions at -^=200 
GeV can well be explained by our model calculation where the continuous regeneration of J /if/ 
from thermalized charm quarks in QGP is an important ingredient. We predict that at higher 
colliding energies, for example at LHC, the regeneration will become the dominant ingredient. 
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